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Abstract: We generalize recent work to construct a map from the conformal Navier 
Stokes equations with holographically determined transport coefficients, in d spacetime 
dimensions, to the set of asymptotically locally AdS^+i long wavelength solutions of Ein- 
stein's equations with a negative cosmological constant, for all d > 2. We find simple 
explicit expressions for the stress tensor (slightly generalizing the recent result by Haack 
and Yarom flarXiv : 0806 . 460^ ) ), the full dual bulk metric and an entropy current of this 
strongly coupled conformal fluid, to second order in the derivative expansion, for arbitrary 
d > 2. We also rewrite the well known exact solutions for rotating black holes in AdS^+i 
space in a manifestly fluid dynamical form, generalizing earlier work in d = 4. To second 
order in the derivative expansion, this metric agrees with our general construction of the 
metric dual to fluid flows. 
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1. Introduction 

The AdS/CFT correspondence establishes a deep connection between quantum field theo- 
ries and theories of gravity. At generic values of parameters both sides of this equivalence 
are complicated quantum theories. However, every known example of the AdS/CFT duality 
admits a large N limit in which the gravitational theory turns classical, and a simultane- 
ous strong 't Hooft coupling limit that suppresses a' corrections to gravitational dynamics. 
In this limit the AdS/CFT correspondence asserts the equivalence between the effectively 
classical large N dynamics of the local single trace operators p n = A _1 TrO n of gauge 
theory and the classical two derivative equations of Einstein gravity interacting with other 
fields. 

The usual rules of the AdS/CFT correspondence establish a one to one map between 
the bulk fields and the single trace field theory operators; for instance, the bulk Einstein 
frame graviton maps to the field theory stress tensor. Given a solution of the bulk equations, 
the evolution of any given trace operator p n {x^) may be read off from the normalizable fall 
off 'at infinity' of the corresponding bulk field. This dictionary allows us to translate the 
local and relatively simple bulk equations into unfamiliar and extremely nonlocal equations 
for the boundary trace operators p n (x). The equations for p n {x) are nonlocal in both space 
and time; indeed the data for the classical evolution of p n includes an infinite number of 
time derivatives of p n on an initial slice. Given the complicated and unfamiliar nature 
of these equations, it is difficult to use our knowledge of bulk dynamics to directly gain 
intuition for boundary trace dynamics. It would clearly be useful to identify a simplifying 
limit in order to train intuition. 

Some simplification of trace dynamics is achieved by focusing on a universal subsector 
of gravitational dynamics |Q] . We focus on two derivative bulk theories of gravity that 
admit AdS^+i X Mi as a solution (here M/ is any internal manifold whose character and 
properties will be irrelevant for the rest of this paper). It is easy to convince oneself that 
every such theory admits a consistent truncation to the Einstein equations with a negative 
cosmological constant. The only fluctuating field under this truncation is the Einstein 
frame graviton; all other bulk fieldsare simply set to their background AdS^+i x Mj values. 
This observation implies the existence of a sector of decoupled and universal dynamics 
of the stress tensor in the corresponding dual field theories. The dynamics is decoupled 
because all p n (x) other than the stress tensor may consistently be set to zero as the stress 
tensor undergoes its dynamics, and this dynamics is universal because the evolution of the 
stress tensor is governed by the same equations of motion in each of these infinite class of 
strongly coupled CFTs. 

While the universal stress tensor dynamics described above is clearly simpler than a 
general evolution of p n {x) in the dual theory, it is still both complicated and nonlocal. It is 
useful to take a further limit; to focus on boundary configurations in which the local stress 
tensor varies on a length scale that is large, at any point, compared to a local equilibration 
length scale (intuitively, 'mean free path') which is set by the 'rest frame' energy density 
at the same point (we will make this more precise below). Local field theory intuition 
suggests that boundary configurations that obey this slow variation condition should be 
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locally thermalized, and consequently well described by the equations of boundary fluid 
dynamics. Hence, we expect the complicated nonlocal T^ u dynamics to reduce to the 
familiar boundary Navier Stokes equations of fluid dynamics in this long wavelength limit. 

All the expectations spelt out above have been demonstrated to be true for d = 4, 5 
from a direct analysis of the Einstein equations (See ||, [| ||[ ||, || ||]). This analysis has 
also been generalized to a large extent for general d in a recent paper by Haack and Yarom 
Q. In this paper we continue and complete the analysis of R in explicating the connection 
between the Einstein equations and fluid equations in arbitrary dimensions . 

In particular, we implement the programme initiated in [Q] to explicitly compute the 
bulk metric dual to an arbitrary fluid flow (accurate to second order in a boundary deriva- 
tive expansion). We verify the expressions for the second order stress tensor dual to these 
flows which was recently derived in @, study the causal structure of the solutions we derive, 
determine their event horizons at second order in the derivative expansion, and determine 
an entropy current for these fluid flows. Further, we compare our results to exact solutions 
for rotating black holes in global AdS^+i and find perfect match to the expected order. 
In the rest of this introduction we will more carefully review some of the closely related 
previous work on this subject in order to place the new results of this paper in its proper 
context 1 . 

The authors of developed a procedure to construct a large class of asymptotically 
AdSs long wavelength solutions to Einstein's equations with a negative cosmological con- 
stant. The solutions in [jjij were worked out order by order in a boundary derivative expan- 
sion, and were parameterized by a four velocity field u^{x^) and a temperature field T{x^). 
These velocity and temperature fields are further constrained to obey the four dimensional 
generalized Navier Stokes equations V M T^ = where the stress tensor T^(x^) is a local 
functional of the velocity and the temperature fields. The form of T^ v (x^) was explicitly 
determined in to second order in a boundary derivative expansion (Some terms in the 
stress tensor were independently determined by the authors of []|]. Especially notable in 
this regard are the pioneering work in 10, pd]].) . Consequently, the construction of Q 



may be thought of as an explicit map from the space of solutions of a distinguished set of 
Navier Stokes equations in d = 4 to the space of long wavelength solutions of asymptotically 
AdSrf+i gravity. 

The spacetimes derived in Q were subsequently generalized and studied in more detail. 
In particular, it was demonstrated in Q that, subject to mild assumptions, these spacetimes 
have regular event horizons. In the same paper, the location of this event horizon in 
the 'radial' direction of AdSs was explicitly determined to second order in the derivative 
expansion and it was found to depend locally on the fluid data at the boundary (via a 
natural boundary to horizon map generated by ingoing null geodesies). The authors of 
also constructed a local fluid dynamical 'entropy current' utilizing the pullback of the 
area form on the horizon onto the boundary. The classic area increase theorem of general 



1 Our main aim here is to provide the appropriate background for our work rather than to review the 
complete expanse of the literature relating hydrodynamics to holography. However, we have included a non- 
exhaustive list of references pertaining to hydrodynamics in the context of holography in the References 
section at the end of this paper. 
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relativity was then used to demonstrate the local form of the second law of thermodynamics 
(i.e., the point wise non negativity of the divergence of this entropy current). On a related 
note, in 0], a formalism was developed for conformal hydrodynamics which describes the 
long wavelength limit of a CFT. Using this manifestly Weyl-covariant formalism, many 
results of and Q could be cast into a simpler form and Weyl covariance could be used 
as a powerful tool in classifying the possible forms of the metric, energy momentum tensor 
and the entropy current. 

In , the construction of [|J was generalized to spacetimes that are only locally asymp- 
totically AdSs, i.e. that asymptote to 

dsl = 1 [dz 2 + d S 2 3>1 ] (1.1) 

where ds 2 1 is an arbitrary slowly varying boundary metric, at small z. It is expected 
that, under the AdS/CFT correspondence, such solutions are a universal subsector of 
the solution space of the relevant CFTs on the Lorentzian base manifold M3 i with the 
metric ds%-±. In agreement with this expectation || demonstrated that long wavelength 



solutions of gravity with asymptotics given by (LI) were parameterized by a velocity and 
temperature field on the manifold M3 1, subject to a covariant form of the Navier Stokes 
equations. As an example of this construction, the authors of § were able to rewrite the 
exact asymptotically global AdSs Kerr black hole solutions in a very simply manifestly 
fluid dynamical form, and demonstrate that the expansion of this metric to second order 
in the derivative expansion is in perfect agreement with the general construction of the 
metrics dual to fluid dynamics at second order 2 . 

All of the results described above were originally worked out for the special case d = 4, 
but some of these results and constructions have since been further generalized. In an early 
paper Van Raamsdonk § generalized the construction of the full second order bulk metric 
to an arbitrary fluid flow on a flat boundary to d = 3 and also computed the holographic 
fluid dynamical stress tensor to second order in boundary derivatives. Some terms in the 
second order stress tensor for the uncharged conformal fluid in arbitrary dimensions were 
calculated using different methods by [|l4|, Further, 1/A and 1/N C corrections to some 
coefficients have been computed in g|, g [Ej| . 



More recently, Haack and Yarom Q partially constructed the second order bulk metric 
to an arbitrary fluid flow in a flat d dimensional boundary (for arbitrary d) and fully 
computed the dual second order fluid dynamical stress tensor for a flat boundary. In 
this paper, we continue the study of Haack and Yarom Q to generalize all of the work 
on solutions of pure gravity duals to arbitrary fluid flows in d = 4 dimensions (reviewed 
above) to arbitrary d > 2. 



2 The authors of |3| also considered the coupling of a slowly varying dilaton to the metric. It would be 
an interesting exercise to consider generalizing the results of this paper to a bulk spacetime with dilaton 
dynamics. However, in this paper, we will confine ourselves to the case where the dilaton is set to its 
background value. On a related note , we should also mention two recent papers^, which appeared 
while this paper was nearing completion in which the fluid gravity correspondence was extended to a class 
of charged black holes in AdSs with flat boundary. 
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This paper is organized as follows. In section 2 below we review the Weyl covariant 
notation for fluid dynamics that we will use in the rest of this paper. In section 3 we 
briefly explain the logic of our construction of long wavelength bulk solutions dual to fluid 
dynamics. In section 4 below we present explicit solutions to Einstein equations to second 
order in the boundary derivative expansion. Our solutions asymptote at small z to 

ds 2 d+1 = ^[dz 2 + ds 2 d _^] (1.2) 

where ds 2 ,_ l 1 is the arbitrarily specified weakly curved metric on the boundary. Our 
solutions are parameterized by a boundary d- velocity field u^{x) and a temperature field 
T(x) where are the boundary coordinates. These velocity and temperature fields are 
constrained to obey the d dimensional Navier Stokes equations, V^T^" = where T^ v is a 
local functional of the velocity and temperature fields. We also present explicit expressions 
for the boundary stress tensor T^ v dual to our solutions. Our answer can be expressed in 
an especially simple and manifestly Weyl-covariant form 



Tfxu = P (9fiu + du^u v ) - 
-2r/r £ 



+ 2rjb 



\t\ , A <7 a/3 cr n . ^ a 



(1.3) 



with 



d 1 
b = — — ; p 



4vrT ' y 167rG A dS& d 

where p is the pressure, T is the temperature, s is the entropy density and r) is the viscosity 
of the fluid. t w denotes a particular second-order transport coefficient of the fluid, a^ u is 
the shear strain rate , lo^ u is the vorticity and C^ap is the Weyl Tensor of the spacetime 
in which the fluid lives. These quantities are described in more detail in section 2 below 
(see also Appendix A for a list of notation used in this paper). Note that our result for the 
stress tensor agrees with those of Haack and YaromQ when restricted to a flat boundary 
manifold, but also includes an additional term proportional to boundary curvature that 
vanishes in flat space. 

In section 5 below , we demonstrate that our solutions all have a regular event horizon, 
and find an expression for the radial location of that event horizon upto second order in 
the derivative expansion. We also construct a boundary entropy current Js that is forced 
by the area increase theorem of general relativity to obey the equation V^Jg > 0. In 
section 6, we rewrite the exactly known rotating black hole solutions in global AdS^+i in 
a manifestly fluid dynamical form. These solutions turn out to be dual to rigid fluid flows 
on S d ~ 1,1 (see pO] for earlier work). These initially complicated looking blackhole metrics 
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admit a rewriting in a rather simple form in the fluid dynamical gauge and variables used 
in this paper. In appropriate co-ordinates, the general AdS-Kerr metric 3 assumes the form 



ds = - 2u^dx^ (dr + r A v dx v ) + r g^ v + u^S v ) X u - oj X u 



dx^dx v 



v^ii ii (1.4) 

+ ' a ^ v dx»dx v 

6 d det [r & - u»„) 

where is the fluid dynamical Weyl-connection and is the Weyl-covariantized 
Schouten tensor introduced in J|]. We demonstrate that the expansion of these solutions 
to second order in the derivative expansion agrees with our general construction of metrics 
dual to fluid dynamics. We end our paper in section 7 with a discussion of our results and 
possible generalizations. 

Note that whereas all the results of this paper pertain apply only to d > 2, these 
results here are easily extended to the d = 2 case which turns out however to be trivial. 
See Appendix B for a brief discussion of the triviality of conformal fluid dynamics in d = 2. 

2. Manifest Weyl Covariance 

In this section we review the Weyl covariant notation we will use in the rest of the paper. 
2.1 Weyl covariant notation for conformal hydrodynamics 

The Fluid-Gravity Correspondence relates hydrodynamic solutions emerging from CFT^ 
to regular Solutions of Gravity in AdSd+i- The defining degree of freedom on the hydro- 
dynamics side are the velocity field u^{x) and the temperature field T(x) and the energy- 
momentum (and an entropy current, see below) are expressed as functions of these basic 
degrees of freedom 4 . The basic equations of hydrodynamics are the conservation of energy- 
momentum supplemented by the constitutive relations which give energy-momentum tensor 
as a function of the basic fields of fluid dynamics. 

In this work, we will in particular focus on hydrodynamics that arises out of a d 
dimensional CFT on a weakly curved non-dynamical metric background gn V . The Weyl 
invariance of the underlying CFT imply that the equations of hydrodynamics are covariant 
under the transformations 

g^ = e 2 ^~g^ , u /t = e^u fl , T = e^f (2.1) 

Quantities (like u^(x)) that transform homogeneously under Weyl transformations are said 
to be Weyl covariant. Notice that ordinary derivatives of Weyl covariant quantities are not 
themselves Weyl covariant. However M developed a 'Weyl covariant derivative', utilizing 



3 Note in particular that for d — 2, cr M „ = = and S M „ term is absent in which case we get the BTZ 
blackhole in AdS3 as shown inQ. 

4 In this paper, we will confine our attention to fluids with no other conserved charge except energy and 
momentum. Fluids with other conserved charges have additional degrees of freedom involving the chemical 
potentials related to those charges. 
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an effective 'Weyl gauge field' built out of derivatives of the fluid dynamical velocity field. 
The formulas of conformal fluid dynamics - and all the spacetime dual metrics to fluid flows 
that we will construct in this paper - may all be written in terms of the Weyl covariant 
derivative of [Q], making their Weyl transformation properties manifest. We will extensively 
use the Weyl covariant notation developed in (2[ in this paper, and so pause to review it 
in detail in this subsection. 

Let us define the 'gauge field' 

A v = u X V x u u - ^-ru u = A V + d v <\>. (2-2) 
a — I 



As indicated by the last part of (2.2), the gauge field A^ transforms like a connection 
under Weyl transformations. This, in turn allows us to define a Weyl-covariant derivative 
on arbitrary tensors of weight w. In particular, the Weyl covariant derivative of an arbitrary 
tensor Qu'.'.'. of weight w (i.e. one that obeys Qu'.'.'. = e~ w ^ 'Qu...) may be defined as 

Qi'r = v A Q';r + w A x Q>;r 

+ [gx*A» - S" x Aa - 5£A X ] QZz + ■■■ (2.3) 
- [g Xu A a - 5^A U - 8% A x ] Q£;; - • • • 

It may be checked that the covariant derivative of a tensor of weight w is also a tensor of 
weight w. 

For example the Weyl-covariant gradient of velocity may be decomposed into symmet- 
ric and antisymmetric pieces as V^Uy = a^ u + uj^ v where is the shear strain rate of the 
fluid and u>^ u is the vorticity tensor. This decomposition demonstrates each of and oj^ v 
are Weyl-covariant with weight w = — 1 (a property that may also be - more cumbersomely 
- directly verified from the definitions of these quantities in terms of ordinary derivatives 
of velocity). 

Recall that in a CFT, the stress tensor transforms with the weight d + 2 upto the 
anomaly corrections and is traceless apart from the conformal anomaly. Using these con- 
ditions, the basic equation of hydrodynamics - the conservation of the stress tensor - may 
be recast into a manifestly Weyl-covariant form 

= V M r^ + ^"(T% - W) = (2.4) 

where W denotes the conformal anomaly of the underlying CFT. 

Below we will have use for additional Weyl covariant curvature tensors that are natu- 
rally constructed out of the Weyl covariant derivative. We will now define these objects 

Pu. , V v \ V x = w Vx + U^x a V a with 
F uu — V fj.Au Wj/Aft 

T^fj,uXa = Rfii/Xa ~l~ F uu9Xo (2-5) 



( A 2 \ 

s Ui9v][*%] ( VaAp + A a Ap - —g a p 
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where -B^] = B IMU — B Ufl indicates antisymmetrisation 5 . As is evident from the expressions 
above, the Weyl-covariantized curvature tensors do not vanish even if the fluid lives in a 
flat spacetime. 

We will also need a few other tensors derived from the Weyl-covariantized Riemann 
tensor including Weyl-covariantized Ricci tensor TZ^ U and Weyl-covariantized Ricci scalar 
T& . In addition, the Weyl-covariantized Schouten tensor S^ u is defined from the Weyl- 
covariantized Ricci tensor in the same way as the ordinary Schouten tensor. 

V = (?V " 5^) = S„ + (V,X + A.A, - f ».,) + ^ (2.6) 

and in turn, the Weyl curvature C^Ao- is related to its Weyl-covariant counterpart C^Ao- 
via 



^fivXa + ^iivgxu 

Conformal hydrodynamics is easily formulated using the tensors above along with Weyl- 
covariant derivatives of velocity and temperature fields. Further, we would also find it 
often convenient to use the following Weyl-covariant combinations 

u X V x a^ = P^P/u x V x a af3 + JZT[^ 

K p „ = + (d - 2) (Vf.Au + A„A U - A 2 g^) + g^VxA x + 

K = R + 2(d- l)V x A x - {d - 2)(d - l)A 2 (2.8) 

V x a^ x = (V A - (d - 1)A X ^ ^a 

V x ujf lX = (V x -{d- Z)A X ) ^ A 

Before ending this section we will now digress to review the motivation for our definition 
of the 'gauge field' in structural mathematical terms Q. The mathematically non inclined 
reader is invited to skip directly to the next section. 

In a spacetime which is associated with a class of Weyl-equivalent metrics, one way to 
write down manifestly Weyl-covariant expressions is to introduce what is known as a Weyl 
connection [21, 23]. A torsionless connection \7 we y l is called a Weyl connection if for 



every metric in a class of Weyl-equivalent metrics, there exists a one form A^ such that 
V™ ey 9v\ = ^Af l g 1 y X . A weyl-connection enables us to define a Weyl-covariant derivative 
acting on a tensor field with a weight 8 w as = V™ eyl + wA^. In terms of £> M , the above 



5 Note that the curvature tensors appearing in this paper are negative of the curvature tensors denned 
in | 

6 defined by the relations TZ^ = TZ l _ l \ L , x and 1Z = 1Z\ X where r R. llv \ a is the Weyl-covariantised Riemann 
Tensor denned in the eqn.(2.5). 

7 We remind the reader that the Schouten tensor is denned by the formula(<i — 2)S' M ^ = i? M „ — 2(d-i) • 
One way to motivate this definition is to look at the definition of the Weyl curvature tensor C^vXa = 

8 by which we mean a tensor field which transforms as Q„ = e~ n "^Q l ^ under Weyl transformation. 
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requirement on the connection becomes T^xg^ = which is just the statement of metric 
compatibility of the Weyl-covariant derivative. Such a mathematical structure is especially 
relevant to the problem at hand since hydrodynamics on a spacetime background provides 
us with a natural Weyl connection. The required A^ is uniquely determined by requiring 
that u^Vxu^ = and T>\u x = 0. Using these equations, we can solve for A^ in a particular 
Weyl frame in terms of the velocity field and the usual 9 covariant derivative V M to get 
the equation ( |2.2[ ). 

2.2 Classification of fluid dynamical Weyl covariant Tensors to second order 

In this subsection, we will classify all the Weyl invariant scalars, transverse vectors (a 
vector /tensor is transverse if it is orthogonal to u^) and symmetric traceless transverse 
tensors with less than or equal to two derivatives. We enumerate these invariants 'on 
shell'. In other words two invariants that are set equal by the equations of motion of fluid 
dynamics are treated as equal in our counting. 

In order to classify the Weyl-covariant tensors, we begin with the basic quantities 
of hydrodynamics - the fluid temperature T = d/(4:irb) and the fluid velocity . The 
temperature is a Weyl-covariant scalar with conformal weight unity and the velocity is a 
Weyl-covariant vector with conformal weight unity. Consequently bu^ is the only derivative 
free conformally invariant object one can build out of the basic fields of fluid dynamics. It 
follows that, at the zero derivative level, there are no non-trivial Weyl-invariant scalars, no 
invariant transverse vector or symmetric traceless transverse tensors. 

Let us now turn to one derivative Weyl invariants of fluid dynamics. Recall that 
we wish to classify only the on-shell invariants. Clearly the relations imposed by the 
VuT^ = on one derivative invariants have their origin in 'zero derivative' or perfect 
fluid contributions to the stress tensor. This tensor is proportional to b~ d (g^v + du^Uy). 
The conservation of this stress tensor allows us to express all derivatives of b in terms of 
derivatives of the velocity(see Appendix C). It follows that, with our rules of counting, at 
the one derivative level, there are no Weyl-invariant scalars or transverse vectors. However 
there is one Weyl-invariant symmetric traceless transverse tensor ba^ u at the first derivative 
level 10 . 

In order to enumerate distinct on-shell Weyl invariants at the two derivative level, it 
is important to account for the contribution of one derivative terms in the stress tensor 
to the equation of motion. As we have explained above, the only possible one derivative 
correction to the stress tensor (in the Landau Lifshitz gauge, see ahead for an explanation) 
is proportional to a^. After analyzing the two derivative relations that follow from the 
conservation of the corrected stress tensor, it is not difficult to show there are three indepen- 
dent Weyl-invariant scalars (which are b 2 a^ u a tl1 ' , b 2 u) llu uJ IJ ' v and b 2 lZ), two Weyl-invariant 
transverse vectors ( bP^T>\a u x and bP^T>\uj v x ) and five Weyl-invariant symmetric traceless 



defined via the Christoffel connection of the metric in that frame 

In this counting, we do not count the pseudotensors since they never occur in our computations. 
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transverse tensors - 

A ~Pfjv ad i A Piw a.B (^-9) 

0> o-\u - j^j (Tapa p and ui^ uj Xu + j^j u: a pu) p . 

We will find this classification of Weyl invariants very useful in the rest of this paper. 
For instance, the classification immediately implies that the energy momentum tensor of a 
general fluid configuration to second order in derivative expansion should assume the form 



T" v = p{g^ + dv?u v ) 



- 2ri - n u x V x a^ + t 2 {uj» x <j Xv + u v \o x ^ 

_1_ C UH P ^ ~<*P~ lit n n.Oinfi ( 2 - 10 ) 

+ U^\^ Xv + -r-u^^ + ... 
a — 1 

where Weyl-covariance demands that 

P OCb~ d , 7/ OC b l - d , 71,2 OC b , 6^ OC b 2 - d (2.11) 

3. Perturbative Construction of Solutions 

In this section, we will briefly review the basic logic that underlies the construction of 
gravity solutions dual to arbitrary fluid flows. The methodology employed in this paper is 
an almost direct generalization of the techniques used in Q [2|, [|, [|, || ||, |^| . Consequently in 
this paper we will describe the logic of our construction and the details of implementation 
only briefly, referring the reader to the references above for more details. 

3.1 Equations of motion and uniform brane solutions 

In this paper we develop a systematic perturbative expansion to solve Einstein's equations 
with a negative cosmological constant 

Gab- d{d ~ 1) G AB = 0, M,N = l...d+l (3.1) 

where Gab denotes the Einstein tensor of the bulk metric Gab ■ 

One solution of these equations is AdS spacetime of unit radius 

ds 2 = — + r 2 (r ]fliy dx fl dx u ) , n,u = l...d (3.2) 

Other well known solutions to these equations include boosted black branes which we write 
here in Schwarzschild like coordinates 

ds2 = TTTT + r2 {-!{r)u il u v dx il dx v + Vnvdx^dx") 

1 d (3 ' 3) 

/(?•) = 9pvU> J, u v = -1, V^v = g^v + u^Uv, h =^f 
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g^ v in (|3.3| ) is an arbitrary constant boundary metric of signature [d— 1, 1), while in the 
same equation is any constant unit normalized d velocity. Of course any constant metric 
of signature (d— 1, 1) can be set to r\n V by an appropriate linear coordinate transformation 
— * hvX v , and can subsequently be set to (1,0 .. . 0) by a boundary Lorentz trans- 
formation. Further b in ( |3.3| ) may also be set to unity by a coordinate change; a uniform 
rescaling of boundary coordinates coupled with a rescaling of r. Thus the d(d + 3)/2 pa- 
rameter set of metrics ([^) are all coordinate equivalent. Nonetheless we will find the 
general form (|3.3| ) useful below; indeed we will find it useful to write ( |3.3|) in an even more 
general coordinate redundant form. Consider 

ds2 = tr^f*^ 2 + f2 (-f(bf)u l ,u v dx^dx u + V^dxPdaA (3.4) 



f 2 /( r ) 



where 



^ = e 2 ^V, u M = e^V 6 = e^6, (3.5) 
(x' 1 ) is an arbitrary function and (L^,, u^, and 6 are as defined in the previous equation. 



This metric is coordinate equivalent to (3.3) under the variable transformation r \— > e~*Y. 
Consequently, the whole function worth of spacetimes ( |3.4|) (taken together with the restric- 
tions ([3.5])) are all exact solutions to Einstein's equations and are all coordinate equivalent. 

While the metrics fl3.4| ) all describe the same bulk geometry, in this paper we will give 
these spacetimes distinct though Weyl equivalent boundary interpretations by regulating 
them inequivalently near the boundary. We will choose to regulate the spacetimes ( |3.4| ) 
on slices of constant f and consequently regard them as states in a conformal field theory 
on a space with metric g^ v {x). With this convention, (the non-anomalous part of) the 
boundary stress tensor dual to ( |3.4| ) is given by 

1 

WnGl dS b d 



T ^ = TTT, Ta (9iw + du^Uu) (3.6) 



which shows that the metric in the Equation ( |3.4[ ) is dual to a conformal fluid with a pressure 
p = l/(167rGAdS^ rf ) an d without any vorticity or shear strain rate. Of course the boundary 
configurations dual to ( |3.4j ) with equal g^ u , u M , b but different values of <fi are related to 
each other by boundary Weyl transformations. 
Notice that 

r >" = F»"' * = -S (3 - 7) 

In other words the velocity field is the unique time like eigenvector of the stress tensor, 
and the inverse temperature field b is simply related to its eigenvalue. We will use this 
observation in the next subsection. 

3.2 Slow variation and bulk tubes and our zero order ansatz 

Consider an arbitrary locally asymptotically AdS^+i solution to Einstein's equations ( j3~l| ) 
whose dual boundary stress tensor everywhere has a unique timelike eigenvector. Let 
this eigenvector (after unit normalization) be denoted by u^(x) and the corresponding 
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eigenvalue by ||. We define to be the d velocity field dual to our solution, and also 

define b(x) to be the inverse temperature field dual to our solution. 

Let 6x(y) denote smallest length scale of variation of the stress tensor of the corre- 
sponding solution at the point y. We say that the solution is 'slowly varying' if everywhere 
Sx(y) S> b(y). (As will be apparent from our final stress tensor below, b(y) may be inter- 
preted as the effective length scale of equilibration of the field theory at y) . Similarly, we say 
that the boundary metric is weakly curved if b(y) 2 R(y) <C 1 (where R(y) is the curvature 
scalar, or more generally an estimate of the largest curvature scale in the problem). 

In the previous section we described uniform brane solutions of Einstein's equations. In 
the appropriate Weyl frame the temperature, velocity, boundary metric and hence the stress 
tensor of those configurations was constant in boundary spacetime. These configurations 
are exact solutions to Einstein's equations. We will now search for solutions to Einstein 
equations with slowly varying (rather than constant) boundary stress tensors on a boundary 
manifold that has a weakly curved (rather than flat) boundary metric. From field theory 
intuition we expect all such boundary configurations to be locally patchwise equilibriated 
(but with varying values of the boundary temperature and velocity fields). This suggests 
that the corresponding bulk solutions should approximately be given by patching together 
tubes of the uniform black brane solutions. We expect these tubes to start along local 
patches on the boundary and then extend into the bulk following an ingoing 'radial' curve. 
However this expectation leaves open an important question: what is the precise shape of 
the radial curves that our tubes follow? 

One guess might be that the tubes follow the lines x^ =constant in the Schwarzschild 
coordinates we have employed so far. According to this guess, the bulk metric dual to 
slowly varying boundary stress tensors and boundary metric is approximately given by 



dg2 = ^ r t^f^ + f2 ( -f{br)u^u v dx^dx v + V^daPda? ) (3.8) 



v ax~ w 

r 2 /( 
where 

9^ = eH 9^, u„ = eV b = e% (3.9) 

where g^uix) is a weakly curved boundary metric, and u^(x) and b(x) are slowly varying 
boundary functions. 

Although this guess seems natural, we believe it is wrong. The technical problem with 



this guess is that the metric of ( p.8[ ) does not in general have a regular future horizon [24] 
(for particular examples of similar metrics that do not have a regular future horizon see [25, 
|26| . [l8|].The last two references show a boost-invariant expansion that develops a singular 
future horizon). In this paper we will be interested only in regular solutions of Einstein 
equations; solutions whose (future) singularities are all shielded from the boundary of AdS 
by regular event horizons. As any perturbation to ( |3.8| ) that turns it into a regular space 
must necessarily be large in the appropriate sense, it follows that ( |3.8[) is not a good starting 
point for a perturbative expansion of the solutions we wish to find. 

There is another more intuitive problem with the proposal that the ansatz ( |3.§| ) is 
dual to boundary fluid dynamics. It is an obvious fact about fluid dynamical evolution 
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that the initial conditions of a fluid may be chosen independent of any 'kick' (forcing) one 
may choose to apply to the fluid at a later time. It seems reasonable to expect the same 
property of the bulk solutions dual to fluid dynamics. 11 Now consider kicking a fluid in 
an arbitrary motion at the point y^. The future evolution of the fluid is affected only in 
the 'fluid causal future' - of y^ . We call this region C(y fJl ). Note that C{y ti ) lies within the 
future boundary light cone of y^ 12 . Now consider the bulk region B(y(fi) that consists of 
the union of all the tubes, referred to above, that originate in the boundary region C(y^). 
Clearly B(y^) is the part of the bulk spacetime that is affected by our kick at y^. Bulk 
causality implies that B(y^) must lie entirely within the future bulk light cone of y^. 

This requirement is not met if we generate B(y^) with our tubes that run along lines 
of constant in Schwarzschild coordinates. However it is met in a particularly natural 
way (given the massless nature of the graviton) if our tubes are chosen to run along ingoing 
null geodesies. 13 

With this discussion in mind, let us consider the ansatz 



2 



ds 2 = -2u^dx^ (dr + r A v dx v ) + r 2 g^ v dx^dx u + -j—yU^Uydx^dx" (3.10) 
where once again 



9^v — ^ ^g^vix 



b = ^b{x), (3.11) 



and g^uix) is a weakly curved boundary metric, and u^x) and b(x) are slowly varying 



boundary functions. When g^ u and b are all constant, ( 3.10 ) is once again simply 
the uniform brane solution, rewritten in Eddington Finkelstein coordinates; i.e. when g^ u 
and b are all constant ( 3.1C| ) and (3.4) are coordinate equivalent (via large co-ordinate 



transformations) . 

However when g^ u and b are functions of x^ ( p.10 ) and ( |3.8| ) are inequivalent and 



in fact differ qualitatively. As we will demonstrate below, under mild assumptions the 
metric in ( |3.10| ) has a regular event horizon that shields all the boundary from all future 
singularities in this space. Consequently, this space may (unlike the spacetime in (|3.4| )) 
legitimately be used as the first term in the perturbative expansion of a regular solution 
of Einstein's equations. Moreover the space described in ( |3.10| ) approximates the uniform 
brane solution along tubes of constant x M in ( |3.10 ); such tubes approximately follow null 
ingoing geodesies in this space. 



For all these reasons, in the rest of this paper we will use (3.1C) as the first term 
in a systematic perturbative expansion of a regular solution to Einstein's equations. The 
perturbative expansion parameter is j4- (we assume that the curvature scale in the metric 
is of the same order as l/5x). We emphasize that the solutions we find could not be 



11 In our set up we can kick our fluid at y M by varying the boundary metric at j/ M (this induces an effective 
force on the fluid). 

12 This is strictly true only if we sum all orders in the fluid expansion. Truncation at any finite order 
could lead to apparent violations of causality over length scales of order 1/T. 

For a related discussion on the desirability of using ingoing null geodesic tubes vis a vis causality 
violating tubes, see |27], Bg] . 
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Figure 1: Penrose diagram of the uniform black brane illustrating the causal Eddington-Finkelstein 
tubes running along ingoing null geodesies . The tubes with ^schwarzschiid~ cons ^ an ^ are a ^ so shown. 
Note that we have suppressed the other regions of the penrose diagram not germane to the discussion 
in this paper. 



obtained in a legitimate perturbation expansion, starting from ( |3.4D . Several authors have 
attempted to obtain the bulk metric dual to a 'boost invariant Bjorken fluid flow' starting 
with the zero order solution described by Janik and Peschanski and correcting it in an 

2 

expansion in l/5xb (that turns into an expansion in 1/ts for those particular solutions). 
As pointed out in p7| , 28|, however, the zeroth order solution of Janik and Peschanski 
is precisely ([O]) for the particular case of boost invariant flow. Consequently, while the 
approach of the current work and ||, ||, ||, [?J are similar in spirit to the perturbation 
procedure initiated by Janik and Peschanski, we differ at a crucial point. While those 
authors effectively adopt (3^) as the starting point of their perturbation theory (for the 
single solution they consider), in our work we adopt the inequivalent and qualitatively 
different space ( |3,10D as the starting point of our perturbative expansion. 



3.3 Perturbation theory at long wavelengths 

The logic behind - and the method of implementation of - this perturbative procedure 
have been described in detail in Q and also in || [|, g, . It has also been described in 
those papers how this perturbative procedure establishes a map between solutions of fluid 
dynamics and regular long wavelength solutions of Einstein gravity with a negative cosmo- 
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logical constant. The discussion in the cited references applies almost without modification 
to the current work, so we describe it only very briefly. 

We start with the ansatz guN = 9mn + e 9MN + ^9mn + • • •• Here gfj N is given by 
( 3.10 ), e is the small parameter of the derivative expansion, and g^ N are the corrections 
to the bulk metric that we will determine with the aid of the bulk Einstein equation. 

In implementing our perturbative procedure we adopt a choice of gauge. As in all the 
metrics described above, we use the coordinates r, for our bulk spaces. We use as 
coordinates that parameterize the boundary and r is a radial coordinate. In order to give 
precise meaning to our coordinates we need to adopt a choice of gauge. In this paper we 
choose the gauge g rr = together with g r ^ = —u^. The geometrical implication of this 
gauge choice was discussed in ||, where it was explained that with this choice lines of 
constant x^ are ingoing null geodesies along each of which r is an affine parameter. Note 
that the gauge choice described in this paper is different in detail from that employed in 
Q and also in |, |, | §. 

The Bulk Einstein equations decompose into 'constraints' on the boundary hydrody- 
namic data and 'dynamical equations' for the bulk metric along the tubes which are solved 
order by order in the derivative expansion. The dynamical equations determine the cor- 
rections that should be added to our initial metric to make it a solution of the Einstein 
equations. At each order, we get inhomogeneous linear equations -but, with the same ho- 
mogeneous parts. These inhomogeneous linear equations obtained from Einstein equations 
can be solved order by order by imposing regularity at the zeroth order future horizon and 
appropriate asymptotic fall off at the boundary. These boundary conditions - together with 
a clear definition of velocity, which fixes the ambiguity of adding zero modes - give a unique 
solution for the metric, as a function of the original boundary velocity and temperature 
profile inputted into the metric g^ N - order by order in the boundary derivative expansion. 

Now, we turn to the 'constraints'. The 'constraints' on the boundary data can be 
shown to be equivalent to the requirement of the conservation of the boundary stress 
tensor. Recall that we have already used the dynamical Einstein equations to determine 
the full bulk metric - and hence the boundary stress tensor - as a function of the input 
velocity and temperature fields. It follows that the constraint Einstein equations reduce 
simply to the equations of fluid dynamics, i.e. the requirement of a conserved stress tensor 
which, in turn, is a given function of temperature and velocity fields. 

It may be worthwhile to reiterate that, as expected from fluid-gravity correspondence, 
metric duals which solve Einstein equations can be constructed only for those fluid config- 
urations which solve the hydrodynamic equations. In the next section, we will present the 
metric which is obtained by adopting this procedure. 



3.4 Weyl Covariance 

In this subsection we explain that the bulk metrics dual to fluid dynamics must transform 
covariantly under boundary 'Weyl' transformations. See || for a more detailed explanation 
of this fact. 
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To start with we note that our bulk gauge choice (described in the previous subsection) 
is Weyl covariant. Any metric that obeys that gauge choice can be put in the form 



ds 2 = -2u ll {x)dx IJ '{dr + V„(r, x)dx u ) + 0^(r, x)dx^dx v (3.12) 

where (5^ is transverse, i.e., u^<5^ v = 0. 14 

For later purposes, we note that the inverse of this bulk metric takes the form 

(3.14) 

+ {&~ l Y v {d^ - v^d r ) ® (a„ - v„$.) 

where the symmetric matrix ((S -1 )^ is uniquely defined by the relations u fl (<5~ 1 ) >11 ' = 
and (<5- 1 y x ® Xu = Si + = 

Consider now a bulk-diffeomorphism of the form r = e _< ^f along with a scaling in the 
temperature of the form b = e^b where we assume that (j) = 4>{x) is a function only of the 
boundary co-ordinates. The metric and the inverse metric components transform as 



V„ = e 



v 



, - e%, = 0^ and ((jr 1 )"" = (d" 1 )"" (3.15) 



Recall however that, within our procedure, the quantities <&av and are each functions 
of and b. Now -u' 4 and b each pick up a factor of under the same diffeomorphism 
(the transformation of b is determined by examining the action of the diffeomorphism 
on (|3.10D ). We conclude that consistency demands that V M and 0^ are functions of b 
and that respectively transform like a connection/remain invariant under boundary 
Weyl transformation. It follows immediately that, for instance <&av is a linear sum of the 
Weyl invariant forms listed in section 2, with coefficients that are arbitrary functions of 
br. Similarly, — rA^ is a linear sum of Weyl-covariant vectors(both transverse and 
non-transverse) with weight unity. 

Symmetry requirements do not constrain the form of these coefficients, which have to 
be determined via direct calculation. In the next section we simply present the results of 
such a calculation. 



4. The bulk metric and boundary stress tensor to second order 
4.1 The metric dual to hydrodynamics 

Using a Weyl-covariant form of the procedure outlined in we find that the final metric 



14 All the Greek indices are raised and lowered using the boundary metric g^ v defined by 

g^u = lim r~ 2 [® M „ - u^V v )} (3.13) 

1 — »oo 

and is the unit time-like velocity field in the boundary, i.e., g^u^Uv = —1. 
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can be written in the form 
ds 2 = -2u^dx^ (dr + r A v dx v ) + r 2 g^ u + u^S^u* - uj^ x uj\ v 



dx^dx v 



+ j^y d {r 2 ~ l -LU a ^)u^u u dx»dx v + 2{br) 2 F{br) 



- 2{br) 2 

- 2(br) 2 H l {br) 



y aj3 



KiM^f^-P^ + Kzibr) 



(br) d 2(d-l) (6r)«* "> ° \ 



dx»dx v 
dx^dx v 



w 



a aj 3(T 



a/3 



dx^dx v 



+ 2{br) 2 H 2 {br) u x V x a^ + a;/<7A„ + ^ A ^ A 



(4.1) 



We have checked using Mathematica that the above metric solves Einstein equations upto 
d = 10. 

The various functions appearing in the metric are defined by the integrals 



i i 

+ 



+ - JL + 

r ft ' 



br d{br) d (d + l)(br) d + 1 (br) 2d 



H 



foo y d-2 _ 1 

1 1 



+ 



1 



+ 



# 



/"Cx 

.{br) = 

Jbr 



- 1 iii" 1;./: 



-F(6r) 2 



2(6r) 2 d(br) d (d + 2)(br) d + 2 (br) 2d 
d£ rt 



+ ... 



y d - 3 dy [1 + (d - T)yF(y) + 2y 2 F'(y)] 
d£ ^ y d ~ 2 - 1 



2(6r) 2 d(fr 
1 



1 roo y d-2 _ i 



+ 



y(y d - 1) 

3d + 5 



d{br) d + l 2{d+l)(d + 2)(br) d + 2 [br) 2d 

Jbr ? Jf 



+ 



2(6r) 2 d(d+l)(6r) d+1 (d + l)(d + 2)(6r) d + 2 (6r) 



2d 



A"; 



(br) ee / 
Jbr 



de 
e 2 



i-^(e-i)F'(e)-2(d-i)e 



+ (2(d - l)£ d - (d - 2)) dy y 2 F'(y) 
(d-3)(d-l) 2(d-2) 1 



2(d + l)(6r) 2 



+ 



# 



d(2d- l)(6r) d (5r) d + 2 



+ ... 
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poo roo _ I 

1 1 

+ 



d(d + 2)(6r) 2 (d+l)(6r) 
1 



r /t \oj ~ • • • 



(d + l)(2d + l)(&r) d+1 2(d + l)(d + 2)(6r) d + 2 (6r) M 



where we have also provided the asymptotic forms of the functions as their arguments tend 
to infinity. 

Later in this paper, we will find it convenient to work with other equivalent forms of 
the above metric. Using 



1 



d-2 



1 



d-2 



x K 

X(J " " 2(d - l)(d - 2)"" + 



(4.2) 



we can write 



ds z 



-2u^dx^ (dr + r A v dx v ) + r 2 g^dx^dx" 



K 



d _ 2 -~- d-2~"~ ' (d-l)(d-2) 



+ tt^j (r 2 - -u^u^Ku^dx* + 2(6r) 2 F(6r) 



2(6r) 2 



+ F(br) 



T a(3 



TS tU ^ C^ a/3 '" M d L£ r/i a apV aP L(br 



dx»dx v 
dx»dx v 
dx^dx h 



2{brf Hiibr) 



a/3 



da^dx" 



+ 2(br) 2 H 2 (br) u x V x o^ + u;/<j Ai/ - ct/kj a , 



dx^dx v 



(4.3) 



or alternatively the metric can be written in the form ( |3,12j ) 



ds 2 = -2 Ull dx^{dr + V v dx v ) + (S^dx^dx" with 



\d-2 M 



2(6r) 



r 2 (l - (br) d ) - l -u a pu«P - {br) 2 K 2 {br) a -^- 



+ ... 



1 



rA » + —2 



V x lo\ - V x a\ + 



n 



2(d- 1) 



-^ a L) x a t 



\d-2 M J 



2(6r) 



r 2 (l - (6r) rf ) - ^ a/3 u^ - (6r) 2 ^ 2 (6r)^^ 
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r P^ u - lo\v 



+ 2(br) 2 F(br) 



-2{br) 2 H l {br) 



+ 2(br) 2 H 2 (br) 



-a^ + F(br)a^a Xh 



T a/3 



2{brfK l (br) <T -^-P, p 
a — 1 



r af3 



d-1 

1 + 



U W 



(4.4) 



Using ( 3.14 ), the inverse metric can be calculated . The tensor (& 1 ) flI/ occurring in 
the inverse metric can be calculated as 



2b 2 

- -^nbr) 

2b 2 

+ A r H 1 {br) 
2b 2 

- -^H 2 (br) 



-a^ - F(br)a^ X (T Xu 



d-l 



T a/3 



u 



a — 1 



(4.5) 



u x Vx^ v + lo»xo av + w V" 



+ ... 



We have checked that results of this subsection agree with the hydrodynamic metric 
duals for d = 4 derived by the authors of |1J] and the d = 3 metric derived in || (in 
order to match our results with older work that was performed in different gauges we 
implemented the necessary gauge transformations). In the next subsection, we proceed to 
derive the stress tensor dual to this metric and compare it against the results available in 
the literature. 

4.2 Energy momentum tensor of fluids with metric duals 

The dual stress tensor corresponding to the metric in the previous subsection is given by 



fpu = V {g^u + du^Uy) - 2-qa^ 



-27JT L 

+ 2rjb 



u X V X (J^ v + lo^oxv + ^u X cr^x 



a a p<j 



af3 



(4.6) 



with 



4vrT 



1 



4vr 16vrG A dS& d " 



and 



1 



16vrG A ds6 d 



oo n ,d— 2 _ i 



y(y° 



dy 



(4.7) 



This result is a generalization to the fluid dynamical stress tensor on an arbitrary curved 
manifold in general dimension d reported in |j| |l[ ^, ||, [l5|, |6| for special values of d and 
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most recently by Q for flat space in arbitrary dimensions. The values of r w for some of 
the lower dimensions is shown 15 in the table |4,2| . 
Comparing this expression with (plOD , we get 

^ = ic = 2t/(ti + r 2 ) = 2r/6 , t 2 = t uj and £ w = (4.8) 

Note that these relations between £0-, £cs r i and r 2 quoted above are universal in the sense 
that they hold true for uncharged fluids in arbitrary dimensions with the gravity duals. 
It would be interesting to check whether these relations between the transport coefficients 
continue to hold against various possible generalizations including the generalization to 
fluids with one or more global conserved charge. 16 

Now, we proceed to compare our results against the results already available in the 
literature. Until now, we have found it convenient to closely follow the parametrisation of 
the stress tensor in M. An alternative parametrisation of the energy-momentum tensor 
was presented in the section 3.1 of || - the parameters r n , Ai 2,3 and k defined there can 
be related to our parameters via the relations 

A 2 

n=T n , r 2 = -— , £ J = 4A 1 , £ c = K (d-2) and £ w = A 3 (4.9) 
zrj 

which in turn gives the value of the transport coefficients as 

Tu = b-T UJ , Ai = — , A 2 = -27/Ta, , A3 = and k = -j— - (4.10) 
which agrees with all the previous results in the literature |L5], |3|, [7| . 

5. Causal structure and the local entropy current 
5.1 The event horizon of our solutions 

Although our assumptions can almost certainly be greatly relaxed, for the purposes of this 
section we specialize to boundary metrics that settle down, at late times to either the flat 

15 More generally, the integral appearing in the expression for r w can be evaluated in terms of the derivative 
of the Gamma function or more directly in terms of 'the harmonic number function' with the fractional 
argument (as was noted in 

b 



d 

IE + ^ L °S T ( Z J 



= — — Harmonic [2/d — 1] 

=2/d 



For large d, t u has an expansion of the form r w /6 = 1/2 — 7r 2 /(3d 2 ) + . . .. 

16 In this context, we would like to note that in the presence of a charge, there are more than one natural 
convention for the definition of the velocity - velocity can be defined as the unit time-like eigenvector of 
the energy-momentum tensor (as we have done in the chargeless case) or can be defined alternatively to 
be the unit time-like vector along the charge current. The former is called the Landau frame velocity and 
the latter is termed the velocity in the Eckart frame. The transport coefficients defined above can depend 
crucially on which of these definitions are used. 

While this work was nearing completion, the authors of [jl2j and jL3| reported independently the transport 
coefficients for a particular class of charged black brane configurations with flat boundaries. Interestingly, 
their coefficients continue to obey £ CT = 2r)(ri +T2) (or equivalently 4Ai + A2 = 277171) in the Landau frame. 
As far as we know, the charge dependence of £c is not known yet. Authors of and jl3| report £ w ^ 
for the charged case in the Landau frame. 
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Value of r u /b for various dimensions 


d 


Value of r w /b - J l v y{yd _ x) dy 


(Numerical) 


3 
4 
5 

6 


U Log3 -573) 
iLog 2 

\ (Log 5 + f y/l - ^ - ^ArcCoth y^) 

i( L °g 3 + 373) 


0.247006. . . 
0.346574. . . 
0.396834. . . 

0.425803. . . 



metric on i?^ -1 ' 1 or the flat metric on S^ 1 x time and to fluid flows that settle down at late 
times to uniform brane configurations on i?^ -1 ' 1 or stationary rotating black holes (studied 
in greater detail ahead) on S x time. See for a discussion on how the dissipative 
nature of fluid dynamics makes the last assumption less restrictive than it naively seems. 

Now the event horizon of our spacetimes is simply the unique null hypersurface that 
tends, at late times, to the known event horizons of the late time limit of our solutions. 
In this subsection we will explain how this clear characterization may be translated into 
an explicit and local mathematical formula for the event horizon within the derivative 
expansion. 

Recall that our bulk metric is written in the gauge g rr = 0, g r ^ = —u^, and conse- 
quently takes the form 

ds 2 = -2u^dx^(dr + V„dz") + &^dx^dx v (5.1) 

where we remind the reader that 05^ is transverse and all the Greek indices are raised using 
the boundary metric g^ u . As we have explained before transforms like a connection and 
0^ is invariant under boundary Weyl transformations. 

Let us suppose that the event horizon is given by the equation S = r — r H (x) = 0. The 
normal vector £4 to this hypersurface is simply the one-form dS = £, A dy A = dr — d^rndx^. 
This one-form - and its dual normal vector - can be written in a manifestly Weyl covariant 
(if slightly complicated) form as follows 

Udy A = dS = (dr + V x dx x ) - fs^ds" 
£ A d A = G AB d A Sd B = - V^d r ) - u^K^d r (5.2) 

where we have introduced two new Weyl-covariant vectors = e~^k^ and = e~~^n^ 
defined via 

^ = d^r H + V^h and 

We use the subscript H to denote that the functions are to be evaluated at the event- 
horizon. 

If we adopt the boundary co-ordinates x^ as the co-ordinates on the event horizon, the 
induced metric on the horizon can be written as 

ds 2 H = [G AB {y)dy A dy B ] r=r {x) = H liV (x)dx^dx v (5.4) 
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with 



AcB 



and the null-condition on the horizon, [Gab]h£, £ = 'H tJi i / n fJ "n L/ = translates to 



(fiS- 1 )""^^ = 2u% 



(5.5) 



(5.6) 



We now follow H to compute the event horizon of our solutions in the derivative 
expansion. We start from a Weyl-covariant derivative expansion for r H given by 



r H = \ + b (h x a afi a aii + h 2 u a poj aP + h 3 K) + . 
= r (o) , r (2) , 

H H ' ■ ■ ■ 



(5.7) 



Note that, since there is no first order Weyl-covariant scalar, 17 there are no corrections to 
r H at the first order in the derivative expansion. 
We first compute 



V^ 1 - S^ x u A - 2L H P»V x a 



+ u. 



K 2H 



a/3 



,(2) 



+ ... 



Substituting the above into (|5.6|) , we get 



r (2) = 26 
H " d 



1 



K 2H a p 



To bring this to the form ( pV7| ) , we use ( p.l|) and ( |4.2| ) to write 



d-d^ )PACT ^ + 

2 



d- 1 



11 • d-2 
2(d-l)(d-2) 



+ 



and 



/OO /"OC 
'd£ y d// 



l 



, c y 3 (y d - 1) 2d 
which gives us the position of the event horizon as 



r H = - + b (hxaapo** + h 2 u a ^ + h 3 1l) +... 



where 



h 2 



2(d 2 + d - 4) K 2H 



d 2 (d-l)(d-2) d(d-l) 
and /13 = 



2d(d - 2 
with K 2 h 



00 de 



d{d- l)(d-2) 

1 -ece- 1)^(0 -2(d-i)£ d - 1 



+2 [id- IX 1 -(d-2)) / dyyV(y) 2 



7 See [q] for a classification of the possible Weyl-covariant tensors. 



(5.8) 



(5.9) 



(5.10) 
(5.11) 



(5.12) 
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5.2 Entropy current as the pullback of Area form 

Once the event-horizon is obtained, one can compute an area form on the horizon which 
when pulled-back to boundary along the ingoing null geodesies gives the entropy cur- 
rent. This general prescription by [|j translates into the following expression for the bound- 
ary entropy current 18 



^GAdS 
^GAdS 



-w 



(5.13) 



(n) , 

where we will define det\\T-l in the following. 

To define detjjjW we will split the boundary co-ordinates x^ 1 to (v, x l ) and we continue 
to use the same co-ordinates also on the event horizon. Under this split, the components 
of the also spilt into (n v ,n l ). We will denote the d — 1 dimensional induced metric on 
the constant v submanifolds of the event horizon by fjjj.Then, we define 



n v \J — detg 



(5.14) 



where is the boundary metric and the expression on the right hand side has been as- 
sumed to be pulled back from the horizon to the boundary via the ingoing null-geodesics. 

(n) 

Though we have used a particular split to define det\\TC, it can be shown that the answer 
that we get in the end is independent of the split(See section 3.3 of Q). Hence, the expres- 
sion in ( |5.13 ) constitutes a specific proposal for what the entropy current of the boundary 
fluid should be. This construction has the advantage that the second law in the boundary 
theory is automatically guaranteed by the area increase theorem in the bulk. 19 

By following the procedure just outlined, the dual entropy current of the conformal 
fluid can be calculated. We get 



Ax a aP a al3 + A 2 u; a(3 uj ap + A 3 K 



a/3 



+ 6 2 



/(A 



(5.15) 



with 



-4-i 



(d + 2) 



K 1H d + K 2H 
d 



An 



1 

2d' 1 



+ ... 



B-2 



d-2 



-2 A? 



(5.16) 



d(d - 2) 



18 A more detailed justification of this formula can be found in Q 

19 The Area increase theorem states that under appropriate assumptions the area of a blackhole can never 
decrease. This statement was proved by Hawking for the case of asymptotically flat spacetimes and is by 
now standard text book material (see e.g. [5^, |3lJ). This theorem has since been extended to black holes in 
more general spacetimes (see e.g. ]32| [33| , |34|]), including asymptotically AdS spaces (see |^| and references 
therein for a clear statement to this effect). 
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where K\nd + K 2 h is given by the integral 



K 1H d + K 2H 



00 d£ 



l-m-l)F\i)-2{d-l)i d - 1 
+2 ((d- l)^ + l) j°° dyy 2 F\yf 



(5.17) 



5.3 Second law and the Rate of entropy production 

In the absence of a clear field theoretic microscopic definition, it may be pragmatic to regard 
the entropy current of fluid dynamics as any local functional of the fluid dynamical variables 
whose divergence is non negative on every solution to the equations of motion of fluid 
dynamics, and which integrates to the thermodynamic notion of entropy in equilibrium. 
According to this characterization the entropy current is any local Boltzmann H function, 
whose monotonic increase characterizes the dissipative irreversibility of fluid flows. 

In the previous subsection we have used the dual bulk description to give a 'natural' 
bulk definition of the entropy current that satisfies all these properties. However, at least 
at the two derivative level, the construction of the previous subsection is not the unique 
construction that satisfies the requirements spelt out in the paragraph above. 

In this subsection we will take a purely algebraic approach to determine the most 
general Weyl covariant two derivative entropy current that has a non negative divergence, 
given the equations of motion derived above. The entropy current of the previous section 
will turn out to be one of a 4-parameter class of solutions to this constraint. 

The most general entropy current consistent with Weyl covariance 20 can be written as 



4G Ad s6 d ~ 1 ^ 



2 a 



b A u 



a/3 



+ b 2 BiV x r 



B 2 V x lo 



(5.18) 



+ 



Since we want to constrain the entropy current upto second order we will need to calculate 
the divergence of this current. In order to perform the calculation in a Weyl covariant 
fashion we note that the ordinary divergence J$ can be replaced by the Weyl-covariant 
divergence T>^Jg with 



V„4 = + (w- d)A ll J p s 



as the conformal weight of any entropy current must be d. 



(5.19) 



Let us now take the divergence of (|5.18| ). We find 
4G Ad S& d ~ 1 2Vs = (d-l)b u^V^ 1 + b 2 u"V^ [Aia afl a af> + A 2 co a ^ + A 3 K 



+ b 2 V 



(5.20) 



20 We assume that there are no pseudo-vector contributions to the entropy current which can possibly 
appear only in the case of d = 4. See M for an analysis in d = 4 including pseudovectors. 
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which can in turn be evaluated using the following identities: 21 



(d - 1) b u^V^b 



-i 



u"X>„ 



a a/3 a 



a/3 



a/3 



pd 

2a^u X V x a'" / 



vPVpTl = -2a^7V + u^F*" + 2V, l V v ^ v - 2{d - 2)V^ [u v F v \ ( 5 - 21 ) 

= -2{d - 2)0^ [a fl x a Xu + u„ x u X v + u x V x a^ v + C^uV 

+ cj^ u + 2V il V v o» v - 2{d - 2)P M [u v F v ] 
(d-3) 



<*V^ 



Substituting for the energy-momentum tensor from ( |4.6| ) and keeping only those terms 
which contain no more than three derivatives, we get 22 



4G AdS fc d - 1 2V£ = ^ 



2A 

<V ~ bd(d - 2) ( A 3 - j-^ ) w m aw\ 



bd(d - 2) (A 3 + d ^_ 2 ^ j (<J^<?\v + « A X>A<V + C mv(3 U a U^ 



(5.22) 



+ (Ai6d + r w ) u a Pa^, 



+ 6 2 (B 1 + 2A 3 )P At P^ + ... 



We rewrite the above expression in a more useful form by isolating the terms that are 
manifestly non- negative (keeping terms containing no more than three derivatives): 



1 [IV 



2b 
~d 

bd(d - 2) 



bd(d - 2) / 2A 2 \ 



d- 2 



1 



d(d - 2) 



+-(A 1 bd + T w )u x V x a^ l 
+ b 2 (B 1 + 2A 3 )V^V u a f * u + . 



(5.23) 



The second law requires that the right hand side of the above equation be positive semi- 
definite at every point on the boundary. This gives us the single constraint : 



#1 + 2,43 = 



(5.24) 



21 The first of these identities is just the re-statement of energy conservation u tl T> v T^ 1 ' = 0. The rest of 
them can be obtained by exploiting the properties of various Weyl-covariant quantities which are detailed 
in Note however that the curvature tensors used here are the negative of those appearing in |J. 

22 Note, in particular, that is zero on-shell upto second order in the derivative expansion(See (C.3)). 
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Equation ( [5.24 ) is the main result of this subsection. Any Weyl covariant entropy cur- 
rent that obeys the constraint spelt out in ( |5,24| ) has a manifestly non negative divergence 
of the entropy current, keeping only terms to the order of interest. 

A particular example of such a current was constructed in [|[. The entropy current 
proposed in pj is equivalent to the following proposal for the coefficients 



Ik 
bd ' 



— , B 2 
2d' 



£1 



-2A ? 



d(d - 2) 



(5.25) 



which yields a simple manifestly non-negative formula for the rate of entropy production 
TV^ = 2 V a^ + .... 

Another example of an entropy current whose divergence is non-negative is the entropy 
current derived from gravity in the previous section using the coefficients appearing in the 
equation ( 5.16 ). Since the values of A3 and B\ appearing in ( |5.16| ) satisfy the constraint 
( 5.24 ) , we conclude that the entropy current constructed in the previous subsection satisfies 
the second law. More explicitly, by substituting the values of A and B coefficients from 
(5.15), we get the rate of entropy production as 



U V J S 



2i] 



1 2 



+ ~{A x bd + t u ) u £>ao> 



+ 



(5.26) 



where Ai and t u have been defined in equations ( 5.15 ) and (4/7) respectively. 



6. Black Holes in AdS 



6.1 AdS Kerr metrics as fluid duals 

In the previous sections, we have found the bulk dual to arbitrary fluid dynamical evolutions 
on the boundary, to second order in the derivative expansion. In this section, we now 
proceed to test our results against a class of exact solutions of Einstein's equations. 

This class of solutions is the set of rotating black holes in the global AdS spaces. The 
dual boundary stress tensor to these solutions varies on the length scale unity (if we choose 
our boundary sphere to be of unit radius). On the other hand the temperature of these black 
holes may be taken to be arbitrarily large. It follows that, in the large temperature limit, 
these black holes are dual to 'slowly varying' field theory configurations that should be well 
described by fluid dynamics. All of these remarks, together with nontrivial evidence for 



this expectation was described in |2C]. In this subsection, we will complete the programme 
initiated in [20| for uncharged blackholes by demonstrating that the full bulk metric of 
these high temperature rotating black holes agrees in detail with the 2nd order bulk metric 
determined by our analysis earlier in this paper. This exercise was already carried out in 
for the special case d = 4. 

Consider the AdS-Kerr BHs in arbitrary dimensions - exact solution for the rotating 



blackholes in general AdS^+i in different coordinates is derived in reference [36]. Following 
|36], we begin by defining two integers n and e via d = 2n + e with e = d mod 2. We 
can then parametrise the d + 1 dimensional AdS Kerr solution by a radial co-ordinate r, a 
time co-ordinate i along with d — 1 = 2n + e — 1 spheroidal co-ordinates on S rf_1 . We will 
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choose these spheroidal co-ordinates to be n + e number of direction cosines fii (obeying 
Efc=i A/e = ^ ) anc ^ n + 6 azimuthal angles with n +i = identically. The angular 
velocities along the different fys are denoted by a, (a n+ iis taken to be zero identically). 

In this 'altered' Boyer-Lindquist co-ordinates, AdS Kerr metric assumes the form (See 
equation (E.3) of the pq| ) 

$dr 2 , 2M f Tirjl ^aitfd<pi\ 2 



ds z 



-W(l + r 2 )di 2 + 



1-2M/V V$ 

n + e 2 _i_ 2 -i 



8=1 



^ r 2 + fl 2 

E Y^f^h 



(6.1) 



where 



^-Er^ ; v^i+^Yia+fy and 

i=l * i=l i=l 

We first perform a co-ordinate transformation of the form 

dr , ajdr 



n+t 0-2 



r 2 _|_ a 2 



dt = dt 



; = dipt 



(1 + r 2 ) (1 - 2M/V) 
followed by another transformation of the form 

-,2 



(r 2 + a 2 ) (1 - 2M/F) 



(6.2) 



(6.3) 



1 



Mi 



to get 



where 



W V 1 - ai 



with 



1 



1 - E< a ?M? 



, 5 = W 



E 



i + 



— rr 1 r" 2 



(6.4) 



ds 2 = -2 Ufl dx»(dr + rA v dx v ) + [r 2 ^ + E^] dx»dx v + ^^dx^dx v (6.5) 



= $ + , A ft = , b = {2M)- 1 / d 

gnu = -eft 2 + ^ (dtf + Hi dip 2 ) 

i 

dt 2 + ^ a? (d/if + i4d<Pi) + ( J] a 2 ^,^ 



(6.6) 



This expression can be further simplified using the following identities 



r 2 V$ = det [r 5% - u" u ] 



where the determinant of a tensor M„ is defined by 



det 



Mi 



(6.7) 
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Hence, we conclude that the AdS Kerr metric in arbitrary dimensions can be rewritten 
in the form 



ds< 



2u fl dx fl (dr + r A v dx v ) + r 2 ^ + u^S^ x u X - u^u^v 



dx»dx v 



+ 



r"u^u v 



6 d det [r 5% - 



7 dx^dx u 



We have checked this form explicitly using Mathematica till d 
This metric can also be written in the form 



ds 2 = -lu^dx^ (dr + r A v dx v ) + ^g^dx^dx" 



1 



(d-l)(d-2) 



dx^dx v 



+ 



or alternatively 



r 2 u^u v 



b d det [r 5% - 



-dx^dx h 



ds z 



-2u^dx^(dr + V u dx u ) + <5^dx^dx v with 



V M = rAfj, - S^u 



9 

r u L 



<5 



2b d det [r 5% - lo^ u ] 



r 2 P^ u - u^uxu 



(6.8) 



(6.9) 



(6.10) 



It is easily checked that this metric agrees (upto second order in boundary derivative 
expansion) with the metric presented in ( [4.3[ ) in section 4 of this paper, upon inserting the 
velocity and temperature fields listed in (|6.6|). 



6.2 The Energy momentum tensor and the Entropy Current for the AdS Kerr 
Black Hole 

The exact energy momentum tensor for the AdS Kerr Black Hole described can be com- 
puted using the standard counterterm methods. The non- anomalous part of the energy 
momentum tensor is given by 



jJLV 



v(9iiv + du^Uy) with p 



l6irG AdS b d 



(6.11) 



which is consistent with (4/7) if we take into account the fact that a^ u = in these 
configurations. 

The equation for the event horizon of the AdS Kerr Black Hole is given by V = 2M or 



i i n 2 

p=<( i +^)II( i +3-) 



r 

H i=\ H 
H 



+ ... 



(6.12) 
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which can be solved for r H to give 



d + 2 



i+£an+... 

b (J^u, 



bU 



(6.13) 



+ ... 



b 2d(d-2)"~ ^ d(d-2)(d-l) 

which agrees with the expression for the event horizon in ( |5.11 ) upon inserting the velocity 
and temperature field configurations (|6.6D . 

The entropy current for the AdS Kerr blackhole can be directly obtained from ( 5.13| ). 
We have the following exact results : 



d<\n = rj" 1 11(1 + 



b d {rl + 1) 



= d t + J2 



r 2 + l 



r l R + a. 



aid,, 



2 ^Z^'ip 



1 - aj o 

r 2 + fl 2 



(6.14) 



±G AdS b d (rl + 1) 



rl + af ^ 



These exact results can alternatively be expanded in a derivative expansion. Keeping 
terms only upto second order in the derivative expansion, we get 



1 r _ d—l, 9 6 2 v^9 1 



TV 



b d ~ 1L ~ d " ' d 
= + b 2 - a 2 ) ai d^ + ... 



(6.15) 



which gives 



4G AdS 6 d - 1 j£ = u" 



2d' 



1 - — u^ UaP 



b 2 lZ 



d(d - 2) 



+ ^P A ^ A (6.16) 



We have checked this form explicitly using Mathematica till d = 8. 

Comparing the above with ( 5.15| ) and remembering that cr a/3 = for the AdS Kerr 
black hole, we find that our results in the previous sections are consistent with these exact 
solutions. 



7. Discussion 



In this paper we have constructed an explicit map from solutions of the (generalized) Navier 
Stokes equations on a d- 1,1 dimensional boundary with an arbitrary weakly curved 
metric g^ u to the space of regular solutions to the Einstein equations with a negative 
cosmological constant that asymptote, at small z, to ds 2 = z~ 2 [dz 2 + g^dx^dx^ . We 
have demonstrated that our solution space is exhaustive locally in solution space. In other 
words consider a particular bulk solution B that is dual, under the map constructed in 
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this paper, to a fluid flow F. Then every regular slowly varying bulk solution to Einstein's 
equations that is infinitesimally separated from B is dual to a fluid flow infinitesimally 
separated from F. 

We have also demonstrated that - subject to certain restrictions on the long time be- 
havior - all the metrics constructed in this paper have regular event horizons, and have 
constructed the event horizon manifold of our solutions in this paper. It would be in- 
teresting to relax the restrictions on long time behavior under which this result follows, 
and simultaneously examine under what conditions these restrictions are dynamically au- 
tomatic from the equations of fluid dynamics. In particular, as the long time limit of a 
fluid flow on a static metric is necessarily non dissipative, it would be interesting to fully 
classify all nondissipative flows on static background geometries. 23 

We have been able to put our construction of the event horizons described above to 
practical use: by pulling the area form on the event horizon back to the boundary, we 
have been able to define an entropy current for the dual fluid flow. The divergence of this 
current is guaranteed to be non negative by the classic area increase theorem of black hole 
physics. The entropy current we have constructed is a sort of local 'Boltzmann H function' 
which can, locally, only be created and never destroyed. The local entropy increase theorem 
establishes the irreversible nature of dual fluid flows. It may be interesting to study the 
structure of gradient flows generated by this 'entropy function'. 

In this paper we have presented explicit expressions, to second order in the derivative 
expansion, for the holographic stress tensor of boundary fluid dynamics. Our expression in 
equation (|4.6| ) , which is only a slight extension of the expression listed in j?]] , is extremely 



simple. It is clear from a glance at equation (4.6) that several ratios of the transport and 
thermodynamical coefficients listed there are 'universal', i.e. independent of the spacetime 
dimension for conformal fluids with a gravitational dual. It would be very interesting to 
investigate whether any of this universality persists for more general (e.g. charged or non 
conformal) fluids with a gravitational dual. It would be useful to have more data in hand 
before speculating further; we hope to return to this issue in the future. 

It would also be useful to understand the second order transport coefficients listed in 
our paper more physically (perhaps in terms of an effective model of microscopic dynam- 
ics that replaces kinetic theory). Along another line, it should not be difficult to derive 
expressions for these coefficients in terms of two and three point correlations functions of 
the stress tensor. In particular expect the expression for the coefficients of the nonlinear 
terms in equation (|4.6| ) to be determined from the three point LSZ-type formula (obtained 
by stripping this correlator of its hydrodynamic poles and extracting the coefficient of k 2 
in the residue). 

It would be interesting to attempt to semiclassically quantize the space of hydrody- 
namic (and related) solutions of gravity with a negative cosmological constant. Such an 
exercise could begin to allow one to account for the effects of statistical fluctuations about 
hydrodynamical flows. 



23 It is natural to guess that this set is exhausted by uniform motion (in the case of the boundary R d ~ 1,1 ) 
and rigid rotations (in the case of the boundary S d_1,1 ), but we are unaware of proofs if any of this intuition. 
We thank G. Gibbons for discussions on this issue. 
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Finally, the map developed in this paper should allow us to incorporate atleast a small 
fraction of the enormous (180 year old) study of hydrodynamics into the study of gravity. 
We hope to return to this issue in the future. 
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Appendices 
A. Notation 

We work in the ( h + . . .) signature. The dimensions of the spacetime in which the 

conformal fluid lives is denoted by d . We usually assume d > 2 unless otherwise specified. 
In the context of AdS/CFT, the dual AdS^+i space has d + 1 spacetime dimensions. The 
greek indices v = 0, 1, . . . , d — 1 are used as boundary space-time indices, whereas the 
latin indices A,B = 0,1, ... ,d are used as the bulk indices. Throughout this paper, we 
take the extra holographic co-ordinate = r with the boundary of the bulk spacetime 
at r — > oo. 

We take the bulk AdS radius to be unity which is equivalent to setting the bulk 
cosmological constant to be AAdS = — d (. d ~ 1 ) _ \y e denote the bulk Newton constant by 
GAdS- F° r t ne ease °f reference, we now give the value of Ga<is f° r some of the well-known 



CFTs with gravity duals : (See f37], |39[ for further details) 



1. The d=4, AA=4 Super Yang-Mills theory on N c D3-Branes with a gauge group SU(N C ) 
and a 't Hooft coupling A = g\ M N c is believed to be dual to IIB string theory on 
AdS 5 xS| =1 with G Ad s 5 = vr/(2iV c 2 ) and ol = (AirXy 1 / 2 . 

2. A d=3, AA=6 Superconformal 24 Chern-Simons theory on N c M2-Branes with a 
gauge group U(N c )fcX U(N c )_fc (where the subscripts denote the Chern-Simons cou- 
plings) and a 't Hooft coupling A = N c /k is conjectured to be dual to M-theory on 
AdS 4 xSj ?=2 /Z fc with G Ad s 4 = A- 2 79A78 = ^-^(IN^/ 2 . 



4 In the case of k = 1, 2, the supersymmetry should get enhanced to d=3, A/"=8. 
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Table of Notation 


Symbol 


Definition 


O 11 

Symbol 


Definition 




Fluid velocity ([T7|) 


T 


Fluid temperature 




Boundary metric 




9nv + u^u v 




b 


d/(4irT) 


P 


Fluid pressure (f4.6[) 




Weyl-Connection (|2.2j) 




Weyl-Covariant 










derivative (|2.3|) 






Weyl-covariantized 




Weyl-covariantized 




T . . /H r-l\ 

Kiemann curvature (^.q) 




Schouten tensor (|2.6|) 


J /II/ 


v^A/] ( see (ELa)) 




See 






Weyl covariantized 




Boundary Weyl 






Ricci tensor/scalar Q2.8|) 




curvature tensor ([2.7D 




Shear strain rate, 1/2 T)r tl u u \ 




Fluid vorticity, 1/2 V^u^ 


V 


Shear Viscosity (|4.q) 




Shear relaxation time ([4,6|) 




Second order transport 


m, K 


Alternate notation for 




Cf ' j /k~\ -i r\l\ 

coemcients (J2.1U|) 


Al,2,3 


2nd order coemcients (|4.9|) 




aiaio n n j • 

AlAdSd+i Bulk metric 


ft 

yAB 


Bulk Einstein tensor 




Defined by $.12) 




Defined by (3.14 


) 




Horizon Normal vector (|5.2|) 




See « 






Induced metric 


r H (x) 


Position of the horizon 




on Horizon (jiOD and (|5.5[) 




given by r = r H ( 


x) 


hi 


Coemcients in (|5.11D 




See $.15) 




F,Hi,Ki,L 


See® 









3. A d=6 superconformal theory on N c M5-Branes is conjectured to be dual to M-theory 
on AdS 7 xS^ =1/2 with G Ad s 7 = 3^ 2 /(16iV 6 3 ). 

These values can be used easily to compare our results with the standard results in the 
literature. 

We use round brackets to denote symmetrisation and square brackets to denote anti- 
symmetrisation. For example, Bt^ = B^ v + B u ^ and -Br^i = B^ u — B vll . Our conventions 
for Christoffel symbols and the curvature tensors are fixed by the relations 

V (t V v = d (t V v + T tl ^V x and [V M , V„]V X = -R^a X V a . (A.l) 

We have included a table with other useful parameters used in the text. In the table [A|, 
the relevant equations are denoted by their respective equation numbers appearing inside 
parentheses. 

B. d=2 

Through the text of this paper we have worked with conformal fluids in d > 2 dimensions. 
In this section we explain that conformal fluid dynamics in d = 2 is special and essentially 
trivial. 
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To start with note that a traceless stress tensor in d dimensions has = d + 2 rf ~ 2 inde- 
pendent components. The assumption of local thermalization in the fluid dynamical limit 
allows us to work instead with the d variables of fluid dynamics; the velocities and tem- 
perature. Now d < s<i for d > 2; it is precisely for this reason that fluid dynamics contains 
physical information beyond the conservation of the stress tensor. However s 2 = 2; conse- 
quently two dimensional conformal fluid dynamics is simply the assertion of conservation 
of the two dimensional stress tensor. One may as well work directly with the components 
of the stress tensor. The general solution to the conservation of the stress tensor in d = 2 
is of course well known. In a frame in which the boundary metric locally takes the form 
ds 2 = e 2< ^dx + dx~ (and ignoring anomaly effects in this discussion) the most general con- 
served and traceless stress tensor is given by T ++ = f(x + ) and = g(x~) for arbitrary 
functions / and g. This constitutes the most general solution to 'conformal fluid dynamics' 
in two dimensions. Note that according to this solution, left and right moving waves do 
not interact with each other. Consequently two dimensional conformal 'fluid' dynamics is 
both trivial and a misnomer; conformal fluids in two dimensions do not locally equilibrate. 

The triviality of conformal fluid dynamics in two dimensions has a simple gravitational 
counterpart: every solution of Einstein's equations in two dimensions is locally AdS^. All 
generally coordinate inequivalent regular solutions of these equations are the BTZ black 
holes. (Note that the point mass solutions, studied extensively for instance in |^(|, have 
a naked singularity atleast from the purely gravitational point of view). Conformally 
inequivalent slicings of the same geometry (a la Brown and Hanneaux) generate the left 
and right moving waves described in the previous subsection. From the bulk point of view 
these solutions are trivial because they are all (large) diffeomorphism equivalent to static 
black holes. 

There is yet another way to express the triviality of conformal fluid dynamics in two 
dimensions. It turns out that there are no non-zero Weyl-covariant quantities which can 
be formed out of velocity/temperature derivatives and hence, as noted by |7|], the first 
order fluid dynamical metric becomes an exact solution of the bulk Einstein equations (see 
section 4 of 0] for more details). For all the reasons spelt out above, in the rest of our 
paper we will focus on d > 2. 



C. Pointwise solution to dynamics at second order in derivatives 

As explained in in order to construct the map from solutions of fluid dynamics to 
solutions of gravity at second order, we need to 'solve' the equations of fluid dynamics, 
at a point x M to second order in derivatives. While it is of course very difficult to find 
the general global solutions to fluid dynamics, the corresponding equations are very easily 
solved at a point. In this Appendix we review the solution of these equations in explicitly 
Weyl covariant terms. The results of this appendix were utilized in our construction of the 
bulk metric in section 4. 

For solving the bulk constraint equations upto second order, we need (T> fJi T tJ,u = 0) 
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evaluated upto second order 



V u b = 2b 2 



4nr] 



T af3 



+ 



(C.l) 



where we have introduced the entropy density s of the conformal fluid related to its pressure 
by s = pd/T = 4irpb. This can be used to solve for the partial derivatives of b completely 
in terms of velocity derivatives 



d u b = A a b + 2b 



2 4tt77 



a/3 



<?a(3Q 

d-1 



+ ... 



(C.2) 



d u d u b = b{d a A v + A^Au) + ... 



Since the left hand side of the last equation is symmetric in ji and u, we get an integrability 
condition 

d„Au = d v A» + ... (C.3) 

Hence, we conclude that to this order we have a valid fluid configuration in a patch 
around a point Pq provided we assume 



b = b Q + eb A u0 x u + 2e 2 b 2 ^- 



a a pa 



d 



H~ 6 ^ -j- *A^,»A/y]Q x^x -\- . . . 



(C.4) 



)1U 



du,A 



+ ... 



For the metric given in the text to be a solution of the Einstein equations, it is necessary 
that the velocity/temperature fields obey the above equations of motion with rj/s = 1/(4tt). 
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